Measurement based entanglement under conditions of extreme photon loss 
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The act of measuring optical emissions from two remote qubits can entangle them. By demanding that a 
photon from each qubit reaches the detectors, one can ensure than no photon was lost. But the failure rate then 
rises quadratically with loss probability. In |T| this resulted in 30 successes per billion attempts. We describe a 
means to exploit the low grade entanglement heralded by the detection of a lone photon: A subsequent perfect 
operation is quickly achieved by consuming this noisy resource. We require only two qubits per node, and can 
tolerate both path length variation and loss asymmetry. The impact of photon loss upon the failure rate is then 
linear; realistic high-loss devices can gain orders of magnitude in performance and thus support QIP. 



Recent experimental successes HHJIEIHI demonstrate that 
measuring two quantum systems can cause them to become 
entangled. Consider an optically active matter system: an 
atom, or an atom-like structure such as a quantum dot or a 
crystal defect. Following laser stimulation such a system may 
emit a photon; its internal state is then correlated with the pho- 
ton. We may generate entanglement by stimulating two sep- 
arate systems simultaneously and monitoring their emissions 
in such a way that we determine characteristics of their mu- 
tual state without learning the source of any given photon [5 J 
(e.g. Fig.[T]l. This approach, and certain comparable alterna- 
tives 1 8 9], have the profound advantage that the component 
systems can be far apart, which may make a scalable technol- 
ogy more attainable. 

Any photon loss from the apparatus will leave the matter 
systems in an uncertain state. Unfortunately real systems will 
always suffer from finite photon capture efficiency, and this 
can be exasperated by non-ideal physics in emission process. 
For example nitrogen-vacancy (NV) defects in diamond emit 
less than three percent of their light cleanly, i.e. without giving 
rise to phonons. A scalable technology would need further 
lossy components for optical routing. Photon loss is therefore 
to be recognised as a fundamental difficulty. 

One solution is weak excitation: reduce the rate of pho- 
ton generation so that emission from both qubits simultane- 
ously is very rare. A detector cUck is taken to imply that 
exactly one photon was emitted, and successfully reached 
the detectors |5 6 7|. However the rate of entanglement 
generation is then inversely linked to the fidelity. Entangle- 
ment between atomic ensembles has been accomplished in 
this way ifTOl . In alternative two-photon schemes both matter 
systems emit a photon, and success is heralded by detecting 
them both UJ] [fl [131 [H [Bl [l6l . We register a failure if ei- 
ther photon is lost, so the success rate falls as the square of the 
transmission probability. Recently two macroscopically sepa- 
rate atoms were entangled in this way [JJ at a success rate of 
about thirty per billion attempts (with a recent preprint report- 
ing a 13 fold improvement |4|). 

Here we describe a new approach to entanglement genera- 
tion that can be fundamentally more rapid. Adopting an en- 
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FIG. 1: Schematic of an apparatus suitable for creating entangle- 
ment between distant entities, using a beam splitter for erasure of 
which path information. Our model incorporates uncertain optical 
path length and asymmetric photon loss via the parameters shown. 
Photon loss is modelled by beam splitters with transmittance Ti and 
T2 ; loss at different parts of the apparatus has an equivalent descrip- 
tion in terms of Ti and T2. We define to be the transmittance 
asymmetry, such that sin(2</>) = (Ti - T2)/(Ti + T2). Asymmetry 
in path length is parametrised by A = ii{x\ — 2:2)/ A. 



tanglement distillation technique, we develop a protocol that 
can quickly implement a perfect gate operation using the im- 
perfect entanglement that results from seeing a lone detector 
'click'. Our gate operation is a parity projection, which sup- 
ports universal quantum computing through the graph state 
approach; we evaluate the performance in this context. The 
approach requires two qubits at each local site. This mod- 
est level of complexity has already been demonstrated in sev- 
eral systems: In atom trap devices one can confine two atoms, 
while recent experiments on NV defects have shown the in- 
teracting electron and nuclear spins are individually control- 
lable ifTTl . The utility of such two-qubit nodes has been stud- 
ied in the context of quantum repeaters 1 1 8 1 . 

In order to provide a clear exposition, we will assume the 
specific energy level structure shown in Fig. [T] Alternatives 
such as A structures can be equally suitable. We begin by 
preparing each of the optically active qubits, which term bro- 
kers following Ref. EOl . in state 



cos(6')|0) +sin(6')|e). 
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Here |e) is the state that decays radiatively to see Fig.[r| 
We collect photons using a lens or a cavity system and di- 
rect them through a beam splitter. Given an ideal apparatus, if 
exactly one photon is detected then the brokers are projected 
onto the (|01) + |10))/ ^/2 state (neglecting any phase that de- 
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pends on which detector clicked, which is trivially corrected 
by a local operation). However in reality we may have non- 
number resolving detectors, photon loss (which may be asym- 
metric) and path length variations. Hence, a single click cor- 
responds to a mixed state 



PB = {1- ^)Z^f |vI/+)(vI/+|Z^-" + r;|ll)(ll|, (2) 

where operator Z'j^'f^ represents the effects of asymmetry in 
the apparatus, as parameterized in Fig. [T] This operator can 
be formally represented as acting on broker qubit Bl alone: 

Zf^f = [cos(0)l +sin(0)ZBi][cos(A)l +isin(A)Z_Bi]. 

(3) 

In Appendix |A] we show that the probability that we indeed 
see a click, i.e. that we obtain ps from initial state |0) is 

PcHck - Tsin'^ie) (2 - rsin2(0) cos2(20)) (4) 

where T — (Ti + T2) /2. The weighting in is shown to be 
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sin^ (61) [2 -Tcos^ (20)] 
2-Tsin2(6')cos2(2(/)) ' 



(5) 



From the work of Bennett et al. II22I we know that there 
is a finite chance of distilling a perfect Bell state from any 
two states of the form (1 - ■q)\^+){^+\ + 77|11)(11|. Note 
this is equivalent to (|2]l without imperfections A and Here 
we show that one can perform a perfect parity projection by 
consuming two (or more) resource states of the general form 
(j2]i. We find that A and can be completely unknown pro- 
vided that they do not drift significantly over the course of the 
protocol, and that two qubits per node suffice for all opera- 
tions. Thus we obtain the high fidelity typical of two photon 
schemes, but at a rate that is affected only linearly by T. 

We present our protocol in terms of a basic iterate that must 
be performed more than once. The aim is to perform a sin- 
gle high fidelity parity projection on two remote client qubits. 
Note that although we are adopting the broker-client termi- 
nology from Ref. |20|, that earlier paper did not exploit the 
brokers for distillation and thus implicitly assumed that the 
initial entanglement mechanism is near perfect (77 and 
Z ^ 1 in Eqn.(|2]l). The first flow chart shows one iterate: 



1. 1 Prepare the two brokers in state \9)\6) and detect 
any emitted photon(s) as in Fig (1). Repeat until 
a photon is detected. Brokers are then in state pg. 



Locally rotate each broker by ( 1 + jtjy )/%/2. 



Perform a local control-phase gate between each 
broker & its own client (a deterministic operation) 



4j Measure both brokers in the X-basis, recording 
"0" or "1" for states |+) and | — ) respectively. 
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Step 2 will map 

|<i>^), where |<f>^' 
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(|00) - |ll))/V2and 



A-^f = [cos(0)l +sin((/))XBi][cos(A)l +isin(A)XBi]. 



Using i and j to denote measurement results on the brokers, 
one complete iterate will transform clients qubits from an ini- 
tial state pc to an (unnormalised) final state (see Appendix [B|: 



Ei.jipc) 



i,3PcPi,jZf-;\ 



-A, 



rjSi^jPcSij, 
(6) 

is a parity projector with parity opposite to 
and Si,j is an undesired projector onto sepa- 



J){i,j\ while P,, 



where Pi j 
that of 

rable states. Specifically Si, 

(1 - {~iy+^ ZciZc2) /2, the labels CI and C2 referring to 
the two client qubits. The variables (pi and A,; are — (— 1)'0 
and — (— 1)*A, respectively. Since this state is mixed, we will 
require least one more iterate. Initially we will consider two 
iterates in total, moving to more general cases later. 



l.| Perform the first iterate, denote measurements 
as ii and ji. 



2. 1 Perform another itCTate, denote measurements 
as 12 and ^'2. 



3. 



Case A 



Compare ii with ^2 and compare ji with j2. 
Case B: Case C: 



Both i and 
J differ. 
Success. 



Either i or j 
differ, not both. 
Failure. 




Neither i or j differ. 
Uncertainty . 



In Case A the we have measurement results of the same 
parity, but differing specific i, j values. The correspond- 
ing parity projectors are consistent, Pi^ j^^Pi^ j^ — Pi^.j^, 
but all other projector combinations vanish, Pi^j^S,, 



'5'i„,i„^ii,ji — Si^^j^Si^j^ = 0. We conclude with cer- 
tainty that the clients qubits have simply been acted on by 
ZpY^'Zpf which is proportional to a pure parity 

projection, as Z^Y^'Z^^''^"^' = cos(20)l. The probability 
of reaching Case A after two iterations is 



Pt«,o = cos2(2(/.)(l-7y)V2. 



(7) 



In Case B the parity of i,j has changed and consequently 
PinJnPiiJi — 0' SO we must conclude that the client qubits 
have been projected into a separable state. This is an unrecov- 
erable failure and the clients must be reset before trying again 
(any prior entanglement with other qubits is lost). In Case C 
the quantum operation on the clients still contains some asym- 
metry operator, {Z^f^'^^')"\ and will also produce a mix- 
ture if all i,j have been the same. We could choose to abort 
here, declaring Case C a failure, which would be a 2 iterates 
only (210) strategy. Later we will show that one can also per- 
sist with further iterations in order to resolve the uncertainty. 

In Fig. |2] we characterise the performance in the simplest 
scenario: the 210 protocol acting on clients initialised to 
|+)|+), thus producing Bell pairs at a rate R. The rate is 
expressed in units of 1/t, where t is the time for a single 
attempt at generating a photon. Rate R is optimised by a spe- 
cific 9 as shown in Fig.|2|a) and detailed in Appendix [P] In 




FIG. 2: Attainable rate of Bell pair production, R. (a) Optimising ex- 
perimental parameter 9 (corresponding to amplitude of | e) in initial 
state 1 6')) to maximise R. (b) Linear and log-log (inset) plots compar- 
ing the rate of Bell pair production in units of 1 /r using our approach 
(green) versus a comparable prior scheme, i.e. the double heralding 
protocol 1141 (blue). Other schemes involving detecting two-photons 
llllll2|[T3lll4| [T51 will have similar performance in the high photon 
loss domain. Note the crossover point of 15.2% above which our 
approach is unhelpful, (c) Regions in which our approach (green) 
and double-heralding (blue) have superior rates of generation, given 
a finite probability of a dark count occurring while we monitor for 
photons following excitation of the matter systems. Here we assume 
a minimum acceptable fidelity of 1 — 10^^, thus there is a 'no-go' 
parameter region where both approaches fail, (d) A top-down view, 
over a wider range of dark counts. For simplicity we assume here 
that photon loss is symmetric: Ti = T2 = T. 



Fig. ^h) we compare our approach to a two photon entangle- 
ment scheme, finding three orders of magnitude improvement 
when T reaches 10~^ (the approximate value that has been 
achieved experimentally [1|). 

Dark counts can be a primary cause of infidelity in en- 
tanglement achieved by the path-erasure approach fT). The 
approach described here, with its single broker-client pair at 
each location, is vulnerable to accumulation of dark count 
noise over successive iterations. However because the num- 
ber of iterations is so small, the accumulation need not be 
substantial, hi Fig. ^c) we make the same comparison as in 
[2|b) but now with finite dark counts and insisting on a mini- 
mum fidelity of 1 — lO^''. The thresholds defining the 'no-go' 
region are nearly identical: the reference scheme can toler- 
ate marginally higher dark counts, about 5% greater, and still 
achieve the target fidelity (note the narrow shelf just visible in 
the 3D plot). Note that to enter the no-go region one could 
introduce further ancilla qubits and combine the protocol de- 
scribed here with schemes like Refs. l|T9l l24l l25l l26l . 
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FIG. 3: The fidelity of the 21 strategy in the presence of drifting 
apparatus cut off at fidelity 1 — 10^^ . Drift in the path length is quan- 
tified byDx = ((xi —0:2) — —3:2))/-^' where unprimed (primed) 
variables represent the initial (final) values. Drift in the photon trans- 
mission probabilities is quantified by Dt = 1 - ^(TpWC^VTl), 
with the same priming notation. 



We have shown that our protocol does not require one to 
determine any difference in path length between the 'arms' 
of the device, nor any asymmetry in the transmission proba- 
bilities; the former will cancel and the later only reduces the 
overall probability of success without degrading the entangle- 
ment. However, this conclusion is only valid if the unknown 
quantities do not drift significantly during a given instance of 
the protocol. If indeed such drift has occurred, but is rela- 
tively slight, then it is straightforward to characterise the level 
of resulting error. In Appendix [C] we perform this calcula- 
tion and find that the infidelity e = {-kD^Y + {Dt/'^.Y + 
Ordcr(Z?j., Dt)*- In Fig. [3] we show the fidelity of the parity 
projection of a 2IO strategy, where drift occures between the 
first and second iterate. Assuming that only path length drift 
is significant, the infidelity will be below lO^'^ provided that 
drift < l/(327r). 

Our protocol has applications far beyond simply generating 
high quality Bell pairs. Because we perform a parity pro- 
jection on the client qubits rather than purifying a specific 
state, we can use the same protocol to subsequently entan- 
gle those qubits with other partners. It is now well under- 
stood that one can thus build up a many-qubit resource such 
as a graph state, and so support efficient universal quantum 
computation IZTl |28l |29l . We will evaluate the performance 
of our protocol in the context of the specific growth strategy 
depicted in Fig. |4] There are now two numbers to optimise, 
the rate of new entanglement generation, and the probability 
that a given entanglement operation will succeed. Thus it is 
interesting to generalise our simple 2 iterates only protocol 
in order that 'Case C need not be abandoned. In fact this 
is straightforward: Over multiple iterates, outright failure oc- 
curs only when the parity of measurement outcomes i, j in the 
latest iterate differs from the previously seen parity. Success 
requires that all outcomes have the same parity and the two 
possible instances of that parity have occurred the same num- 
ber of times. The latter criterion assures that any errors de- 
scribed by Z will not degrade the resulting entanglement; our 
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FIG. 4: A suitable method for full scale quantum computation us- 
ing the well studied graph state entanglement. The approach de- 
picted here provides the context for us to evaluate our protocol. 
Each dot (or 'node') represents a qubit, each line (or 'edge') rep- 
resents entanglement between qubits. On the left side the graph state 
has a regular ordering, using which one can perform a computation 
purely by measuring out individual qubits (each in an appropriate ba- 
sis) 1 27, 28, 29 1 . The graph is thus consumed from the left, driving 
the logical qubits to the right such that they are forever on the left- 
most fringe. Simultaneously a second set of measurements, i.e. those 
corresponding to the entanglement process described here, generate 
fresh graph state structure (right edge) and fuse it into a regular ar- 
ray (central). Entanglement is created 'just in time' for use. As each 
physical system is measured out on the left, it is reused in the pro- 
cess of new growth. The small figures to the right show the results 
of success (creation of new entanglement) or failure (destruction of 
prior entanglement). 



former expression for the state projector simply generalises to 
(Z^Y'^')?(Z^^"^'^*)ti-'ij which reduces to a pure par- 
ity projection as before. Given a system where Z ^ 1, this 
criterion could of course be relaxed. 

In the paradigm illustrated in Fig. [4] one would employ the 
simple 210 strategy for the creation of the Bell pair reser- 
voir, and the more sophisticated multi-iterate protocol for the 
creation of other 'edges' where failure will cause damage (as 
shown in the Fig. [4] insets). The characteristic time required 
to perform operations on the logical qubits, i.e. the effective 
clock speed of the quantum computer, depends on the time 
needed to create new edges in the graph. It is straightforward 
to determine this variable in the limit of low photon capture 
probability T; the calculation appears in AppendixjE] The av- 
erage time required to increase a branch length by one qubit, 
accounting for both successes and failures of the form shown 
in Fig. |4ta), is found to be 29.0(r/T). This is in contrast to 
Order(r/T^) for two-photon schemes such as Ref. lfT4l . We 
emphasise that although this 'clock speed' depends on T, in 
our approach the probability of success is substantially de- 
coupled from T. Finally we note that although this example is 
2D, one can equivalently generate higher dimensional graphs. 



including the structures suggested by Raussendorf et al. for 
fault tolerant computation with a high threshold |30|. 

We have described a form of measurement-based entangle- 
ment generation that can be used in situations of high photon 
loss. Performance can be orders of magnitude higher than cur- 
rent experiments, thus bringing fault tolerance thresholds cor- 
respondingly closer and enhancing the prospects for a highly 
scalable technology. We thank Joe Fitzsimons, Sean BaiTett, 
Jason Smith, Pieter Kok and David Moehring for helpful com- 
ments. This research was supported by the Royal Society and 
the QIP IRC. 



APPENDIX A: A DERIVATION OF 77 & PcUck- 

We denote PcUck as the probability of single detector click 
from two brokers in the state \0)\6). This is simply the sum of: 
the probability of a single photon being emitted and captured, 

P(i) = 2Tsm^{e) cos^{e); (Al) 

and the probability of two photons being emitted and either 
being captured 

P(2) =sin4(0)[l-(l-T)2]; (A2) 

Note that the either is inclusive since photons bunch and de- 
tector are assumed to be non-number resolving. The result- 
ing mixed state has a |11)(11| contribution with a proportion 

V = P{2)/Pclick- 

APPENDIX B: THE QUANTUM OPERATION (pc). 

We derive the quantum operation for measurement results i 
& j by considering how each pure state |^'+) and |11) affects 
the clients, and then weighting these outcomes appropriately. 
Considering Z^j^|*+) first, the local rotation (1 + iY)/y/2 
maps this to the state Xgf'~'^\^~). Next a control-Z is per- 
formed between broker and client, followed by an X— basis 
measurement, resulting in: 

(±\bCZ§ = ((0|b ® 11c ± ® Zc)/V2. (Bl) 

It is easy to verify that this operator maps the asymmetry op- 
erators between brokers and clients, such that: 

{MsiCZ^lX^t^-^ = {±\BiCZ^lZ^t'^^. (B2) 

Using iSzj to represent the measurement outcomes on brokers 
Bl and B2, we conclude that the resulting effect on the client 
qubits is: 

{{-iy,{-iy\Bi.B2cz§lcz§lx^f'-^\<^>-)(B3) 
= {{-iy,{-iy\Bi.B2cz^icz§i\^-)ztt' 

= (lcilc2 - i-iy+' ZciZc2)2t'{^^ /2V2 

— ^ij^ci /^^ 
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where (pi = —(—!)*(/) and similarly = — (— 1)*A. Hence, 
the 1^'+) component projects the clients into a parity subspace 
that is the opposite of the measurement outcome, which is a 
projector that we denote Pij. 

When we consider the |11) component, its separability 
means that it is sufficient to derive the effect of a single bro- 
ker on its clients, and the same will hold for the second broker. 
Forasingle brokerthe rotation (]l+iF)/V2 maps |1) |+), 
which is followed by a control-Z operation to give: 



m + \i)Zc)/V2 



(B4) 



so measuring a broker in the state )) means that its 

client is projected into the separable state |0)(|1)). Taking 
both brokers in consideration, the measurement signature i, j, 
results in the projector Sij = \i, j) {i, j\. 

To construct a quantum operation we begin by considering 
that 1 5'+) and |11) contributions will occur with probability 
(1 — 77) and T]. However, the occurrence of a parity projector is 
shared between two different measurement signatures, which 
introduces a factor of 1/2. This gives an (unnormalised) quan- 
tum operation of: 



(B5) 



Note that the probability of each outcome can be calculated 
by taking the trace of Eij (pc)- As a check it is reassuring to 
calculate Ei,je{o,i}tr[i;i,j(pc)] = 1- 



APPENDIX C: DRIFT EFFECTS 

Since ^'?i.^^-<^ oc 1, we know that for two distillation 
iterates with different parity outcomes asymmetry operators 
vanish. Of course, this assumes that the asymmetry is constant 
over the time between two iterates. In this section we consider 
the errors caused by drift or jitter in the apparatus: 



1 



|(*+|*drift) 



drift I Adrift 



t)' 



where. 



I Adrift) = Z't'+^'i>'^+s^Z- 



(CI) 



(C2) 



where Scf) and 6 A quantify the apparatus drift between first 
and second iterate. This expands out to a function of asym- 
metry variables ^ & A, and their associated drifts S(f) & 5 A: 



cos2(20 + d(j)) sin^idA) + sin^iScp) cos2((5A) 
cos'^ {2(1) + 5(j)) + sm^ {S(t)) 



(C3) 



This equation does not directly lend itself to an interpretation 
in terms of physical properties of the apparatus, so we intro- 
duce new variables: 



Dt 



{{X-L-X2) - {x[ -X'2))/X 



(C4) 



= 1- 



It{T2 



where unprimed (primed) physical parameters describe the 
apparatus for the first (second) iterate of distillation. In terms 
of these variables, we find: 



e = 



(C5) 



Expanding this expression to leading order in and Dt, 
we have a small-error approximation that is quadratic in both 
variables: 

e = {7rD^f + {DT/2f + 0[D^,DT]'' (C6) 



APPENDIX D: RATE OF EPR PRODUCTION 

In this section we detail how to calculate the rate that Bell 
pairs can be produced. In calculating this rate we have the 
preparation angle ^ as a free parameter that we may optimise 
over. This varies as a function of the overall photon loss, and 
also the asymmetry in photon loss. 

Using the 210 strategy, the rate of Bell pair production is: 



R 



1 



-ftwo * Pclick/'^-) 



(DI) 



where Ptwo is the probability of succeeding after 2 iterations, 
and Pciick is probability of succeeding at making the noisy 
Bell pairs required to drive the distillation. The factor of 1 /2 
is present since 2 iterations are required for each attempt. Fi- 
nally, the variable r is the time taken for each attempt at the 
optical protocol used to make noisy Bell pairs (i.e. each at- 
tempt within step (I) of the iterate). We neglect the time 
taken to perform the local operations used in the distillation 
protocol, since we are principally interested in the regime of 
heavy photon loss where the dominant time cost comes from 
the many attempts at obtaining ps- Expressions for the prob- 
abilities, Ptwo and Pciick, can be found in the main paper, and 
using these we can obtain the production rate: 



R 



Tcos^(2(/)) sin^(6') 008^(6*) 
2-Tsm^{9) cos2(2(/)) 



/r. 



(D2) 



It is easy to show that in the limit of vanishing Tcos^(2(/)), 
the optimal of our free parameter is sin^(^) = 1/3. For com- 
parison, the strong excitation used in comparable two photon 
schemes would correspond to sin^(6') — 1/2. Since we are 
primarily interested robustness against photon loss, the figures 
provided in the main paper are optimal assuming that there is 
no asymmetry (j) = and that the amount of photon loss is 
known. If these quantities are not known then the effect on 
rate will only be slight. 



APPENDIX E: RATE OF CHAIN GROWTH 

In this section we show how to reproduce the optimal 
growth rate when using a loop strategy to add Bell pairs to 
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a chain. The essential difference between Bell pair growth 
and chain growth, is that the latter requires consideration of 
how distillation failures reduce chain length. This considera- 
tion means that a loop strategy gives a superior efficiency, as it 
avoids failure at any cost. The rate at which qubits are added 
to a chain is: 



G 



(El) 



where the new variable Pioop and (/), respectively represent: 
the probability of success without having to reset the client 
qubits; and the expected number of iterates until a distillation 
attempt is concluded by a success or reset. The two terms in 
the numerator of G are the success and failure contributions to 
chain length. Since we are primarily concerned with robust- 
ness against photon loss, we will neglect the effects of photon 
loss asymmetry. 

Before calculating the success probability, Pioop, and ex- 
pected number of iterates until an attempt is concluded, (/), 
we must first find the various probabilities for every fc*'* iter- 
ate. Denoting k*^ iterate failure and success probabilities by 
Pf(k) and Ps{k), respectively, it follows that: 

oo 

Pioop = ^-P.(/^), (E2) 

fc=2 
oo 

(/) = ^fc(p,(fc) + py(fc)). 

fc=2 

Both success and failure on a k*^ iterate, may occur via a 
number of different sequences of measurement results, which 
must be summed over, so: 



Pfik) 



(E3) 



2Nf{k)ri 



1 — rj 



fe-i 



and. 



Ps{k) = 2N,ik) 



(E4) 



The failure outcome has two terms representing: the probabil- 
ity of faiUng after a sequence of prior iterates with identical 



measurement results; and the second represents the all other 
possible combinations of prior measurement results, of which 
there are Nf{k). These terms differ because in the first in- 
stance, the density matrix is mixed prior to failure, whereas in 
the latter case the density matrix is pure prior to failure. As for 
the probability of success, this again depends on the number 
of different sequences of measurement results, Ng{k), which 
will obviously only be non-zero for even values of k. 

Lastly, we must find a method of calculating the differ- 
ent valid combinations of measurement results, Ns{k) and 
Nf{k). We begin by defining a vector v{k) which describes 
number of different fc — 1 measurement sequences prior to 
failure. Prior to failure the all measurement results will be of 
the same parity, but this leaves two possible outcomes. Us- 
ing the variable i to denote the difference in the number of 
the two possible outcomes, the (i + 1)*'* entry of v(fc) is the 
corresponding number of possible sequences with i different 
measurements. 

Having defined v(fc), we calculate it iteratively, by using: 



v(A:) = M''-\{2), 



(E5) 



where M has non-zero elements everywhere except above 
and below the diagonal, where the elements are 1; Formally, 
Mi,j = 5i,j±i. The vector v(2), at which the first success can 
occur, has only the first entry as unity and the rest are zero. In 
matrix form, we have: 



v(fc) 



/O 1 

' 1 1 

1 

1 









Vo ... 1 0/ 



fc-2 



















(E6) 



Now Ns{k) and Nf{k) follow directly, as the former is sim- 
ply the first component, Ns{k) — vi(fc), and the failure vari- 
able is a sum from i = 2 to the penultimate non-zeo entry. 

We have outlined all the necessary steps required to calcu- 
late G, and doing this for different values of rj will allow use 
to find the optimal growth rate. Performing this calculation 
for small photon capture, T, and negligible capture asymme- 
try, one finds that the optimal G is G = 0.0345r/r. The 
reciprocal of this rate, 29.0r/T, gives the average time taken 
to extend the chain by one qubit. 
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